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ABSTRACT 

Ca ii /i/13934, 3969 absorbers, which are likely to be a subset of damped Lyman a systems, are 
the most dusty quasar absorbers known, with an order of magnitude more extinction in E(B-V) 
than other absorption systems. There is also evidence that Can absorbers trace galaxies with 
more ongoing star-formation than the average quasar absorber. Despite this, relatively little 
is known in detail about these unusual absorption systems. Here we present the first high 
resolution spectroscopic study of 19 Can quasar absorbers, in the range 0.6 ^ z a b s ^ 1.2, 
with Ca ii /13934 equivalent widths, W ( 3934 ^ 0.2 A. Their general elemental depletion patterns 
are found to be similar to measurements in the warm halo phase of the Milky Way (MW) and 
Magellanic Clouds interstellar medium. Dust depletions and ff-enrichments profiles of sub- 
samples of 7 and 3 absorbers, respectively, are measured using a combination of Voigt profile 
fitting and apparent optical depth techniques. Deviations in [Cr/Zn] ~ 0.3 + 0.1 dex and in 
[Si/Fe] > 0.8 + 0.1 dex are detected across the profile of one absorber, which we attribute 
to differential dust depletion. The remaining absorbers have < 0.3 dex (3<x limit) variation in 
[Cr/Zn], much like the general DLA population, though the dustiest Ca n absorbers, those with 
Wq 934 > 0.7 A, remain relatively unprobed in our sample. A limit on electron densities in Ca n 
absorbers, n e < 0.1 cirT 3 , is derived using the ratio of neutral and singly ionised species and 
assuming a MW-like radiation field. These electron densities may imply hydrogen densities 
sufficient for the presence of molecular hydrogen in the absorbers. The Ca n absorber sample 
comprises a wide range of velocity widths, AV90 = 50 - 470kms~ 1 , and velocity structures, 
thus a range of physical models for their origin, from simple discs to galactic outflows and 
mergers, would be required to explain the observations. 
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1 INTRODUCTION 

Quasar absorption line systems are understood to be probes of 
galaxies not subject to the luminosity bias of traditional surveys, 
being instead biased by their gas content and galaxy cross-section. 
Despite this, our exact understanding of how absorbers and galax- 
ies are related is uncertain due to the difficulty of interpreting ab- 
sorption line characteristics in the context of the associated galaxy. 
Damped Lyman a systems (DLAs), defined by virtue of having 
neutral hydrogen column densities Whi ^ 2 x 10 20 atoms cirT 2 , 
are some of the best studied absorption systems. In particular, it 
was postulated that DLAs may be the reservoirs of neutral gas 
in the Universe from which stars form because the amount of 
neutral gas in DLAs a t z ~ 3 is similar to the amount of gas 
locke d up in stars today JWolfell993klWolfe. Gawiser & Proch aska 
2005). Attempts to measure the star-formation associated with 
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DLAs through a variety of direct and indirect techniques (e.g. 
Wolfe. Prochaska & Gawiserl 120031 : IChen. Kennicutt Jr. & Rauch 



2005; Wolfe & Chen 2006) have met with limited success, gener- 



ally finding low star-formation rates. 

Recently. Iwiid & Hewettl d2005h used the statistical power of 
the Sloan Digital Sky Survey (SDSS) to study the rare, so-called 
strong Can abso rbers (those with rest- frame equivalent widths, 
W 3934 > 0.2 A). [Wild & Hewettl 120051) postulated that Can ab- 
sorbers were an unusual sub-class of DLAs because lines of sight 
through our own Galaxy with similar strengths of Can absorbers 
have neutral hydroge n column den sities consistent with DLAs. 
Furthermore, Iwild & Hewettl ( 120051) found that the Mg 11 equiv- 
alent widths of Can absorbers are large (W^ 7% > 1.0A), thus 
they have a high probability of being DLAs teao & Turnshekl 
2000). The observed sample of Ca 11 absorbers has also been found 
to be much more dusty than either the general DLA popula- 
tion or strong Mgn absorbers , with avera ge E(B- V) ~ 0.1 mag 
dWild. Hewett & Pettinil [20061 henceforth |WHP06|) . compared to 
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E(B- V) 0.013 mag for Mgn absorbe rs dWHPOd lYork et all 
20061) and E(B-V) < 0.02 mag for DLAs (Murphy & Liskell2004 



< 0.5 A, thus this can be used as a constraint to weed out 



Vladi lo et a D l2008h . The low levels of dust in DLAs are con- 



firmed by studies of radio-selected quasars, which find little evi- 
dence for optical quasar drop-ou ts due to dusty DLAs ( Ellis on et al.l 
l200lMEllison. Hall & Lirj2 005). Nonetheless, the presence of dust 
detected by IWHPOrj means that Can absorbers are significantly 
under-represented i n existing samples of DLAs, a conclusion con- 
sistent with those o flNestor et al J 12008). who find that only ~10 per 
cent of DLAs possess very strong, W^ 934 ^ 0.5 A, Can absorption. 
The suggested number of optical quasar drop-outs due to dusty Ca n 
absorb ers is still consistent with the statistics of t he lEllison et all 
j200ll) radio-selected quasar sample. iNestor et al.l j2008h also find 
indirect indications from photo-ionisation modeling, constrained 
by their observations, that Can absorbers host molecular hydrogen. 

In addition to these indirect indicators of star-formation ac- 
tivity, such as dust and molecular hydrogen, direct studies at both 
intermediate and low redshift have shown Can absorbers to be as- 
sociated more strongly with star formation than other absorbers. 
At z ~ 1 Can absorbers are found to have mean in situ star for- 
mation rates (SFRs) from 0. 11 - 0.48 M Q yr -1 compared to just 
0. 1 1 -0. 14 M yr 1 for DLAs dWild. Hewett & Pettinil2007l) . These 
SFRs are as measured inside the 3" diameter SDSS fibre centred on 
each quasar and thus do not account for star-formation in the ab- 
sorber host galaxy if it falls outside this radius. At lower redshifts, 
z < 0.4, direct measurements from the host galaxies of Can ab- 
sorbers give SFRs in the range 0.3 - 30 M Q yr~', t hough it is likely 
the work is biased towards L* galaxies dZvch et alj2007l) . 

Thus, there is now much evidence to suggest Can absorbers 
are associated with star-formation, at least to some extent, both with 
direct measurements of the star formation rate and through sig- 
natures generally associated with star formation such as dust and 
molecular hydrogen. However, few details of why Can absorbers 
might trace star-formation more readily than other absorbers are 
understood. Nor do we know much about the physical environ- 
ment of the absorbers themselves (rather than their host galaxies). 
In this paper we present the first results of a high resolution study of 
the absorbing gas using observations from the Ultra- Violet Echelle 
Spectrograph (UVES) at the VLT and the High Resolution Echelle 
Spectrograph (HIRES) at Keck. With these data we are able to 
study the dust depletion and chemical uniformity (a enhancement) 
of individual absorbers as well as constrain the physical environ- 
ment giving rise to the absorbing gas based on the velocity structure 
of the absorber. 

The paper is structured as follows: in Section[2]we present the 
sample selection, observations, data reduction and analysis tech- 
niques used in the rest of the paper. Section [3] contains the results 
of our analysis, divided into dust depletion, enrichment and envi- 
ronment sub-sections. Section [4] places these results into context 
with more general results from the literature and finally, Section|5] 
summarizes the main points from the paper. 



false-positive detections dWild & Hewettll2005l ; Iwild et alj|2006h . 
Absorbers were identified by searching for features with rest frame 
equivalent width, ^ 0.2 A at 4 a significance in the sample 

of all 58835 SDSS DR4 quasars classified by having 'specClass' 
= 3 or 4 in the SDSS. To confirm a detection as real, it was also 
required that the Caii/13969 line equivalent width had a signifi- 
cance > 1 a and that Mg n /12796 was detected at > 6 cr significance 
with W™ 6 ^ 0.5 A. This requirement implies a minimum velocity 
width ~ 50 km s -1 for our absorbers, which being at the lower limit 
of the velocity width distribution for DLAs will not strongly bias 
any comparison between Can absorbers and DLAs (See Section 
14.4b . Absorbers separated by less than 500kms~' were classified 
as a single absorption system. For the purposes of this project the 
sample was restricted to z abs > 0.73 so that Znn/12026 was above 
3500 A, below which the efficiency of UVES (and HIRES) begins 
to drop more rapidly. Finally, suitable candidates for follow-up high 
resolution observations required the Zn and Cr transitions to lie out- 
side the Lya forest. 

Targets with g-band magnitudes ^ 18.5 mag were se- 
lected to achieve signal-to-noise ratio (SNR) ~ lOperpix in 
the continuum at the wavelength of Zn n ,1/12026, 2062 and 
Crn/12056, 2062, 2066. Declination was restricted to allow obser- 
vations from the VLT at Paranal in Chile (except for the one target 
observed with Keck in Hawaii) and RA. was scheduling depen- 
dant. 

These SDSS-selected absorbers were complemented by a 
search through the ESO/UVES science archive for any absorbers 
which also exhibited strong Can absorption at the limits defined 
above. 



2.2 Observations 

The new UVES data were taken in service mode during the period 
2007 May 18 - July 25 as part of ESO observing programme 79. A- 
0656. The programme was not completed, so we only have partial 
data on some absorbers, therefore we will highlight in each section 
which absorbers have data which contributed to the respective anal- 
yses. We observed J0334-07 1 1 on 2008 January 12 for one 2700 s 
exposure with the CI decker on Keck/HIRESb dVogt et aljfl994r) . 
The observations are summarised in TableQ] The emission redshifts 
are as given by SDSS and the absorber redshifts are as measured 
from the SDSS Can detections. For all UVES observations the slit 
width was 1.2", which was comparable to the seeing in all cases. 
The data were binned, 2 x 2, to improve the SNR. All observations 
were taken with the slit at the parallactic angle to compensate for 
atmospheric dispersion. 

To supplement these newly observed data we include 12 Can 
absorbers detected in quasar spectra from the ESO/UVES sci- 
ence archive (See Table |2j. Equivalent width measurements of 
Can/13934 and Mgn ,12796 for all the absorbers are available in 
Table E] 



2 SAMPLE SELECTION, OBSERVATIONS, DATA 
REDUCTION AND ANALYSIS 

2.1 Selection 

The Caii/L13934, 3969 absorption lines are relatively weak and 
thus detection of the Mg n ,1,12796, 2803 absorption lines at the 
same redshift are required to confirm that the absorption sys- 
tem is real. In practice no real Can absorbers are detected with 



2.3 Data reduction 

The UVES data were reduced using the new UVES Common 
Pipeline Language (CPL) pipeline from ESO. The Keck/HIRESb 
data were reduced using the HIRedux software packag^] bundled 

1 http://www.ucolick.org/~xavier/HIRedux/index.html 
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Table 1. A summary of the observations taken with UVES during ESO Period 79 and HIRES for this work. We define the resolving power, R = AA/A, where 
A is the UVES dichroic central wavelength measured in ran. 
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1 HIRESb target which used the CI decker with FWHM ~7 km s 1 spectral resolution 



Table 2. A summary of spectra taken from the ESO/UVES Science Archive for this work, with an indication of how the data is used. Magnitudes are SDSS 
g-band values unless otherwise referenced. B1950 names are included for convenience. 



Object a(J2000) <5 (J2000) B1950 magnitude z C m Zabs Notes' 



J0004-4157 
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Letters denote the fo llowing: (V)elocity structure an al ysis, (D)epletion pr o file analysis, (A)lpha e nhancement analysis 
2 lRao. Turnshek & Nested UwSt 3 lMaza etalU 19931) 4 lReimers et al] 1 19981) 



within the XIDL software packag^E The algorithms calibrate, sky 
subtract and optimally extract the ec helle orders of the spectrum 
(Ber nstein. Buries & Prochaskall2008ri . 

The individual orders and exposures were then air-vacuum 
and heliocentric corrected and combined using uves popler0 Dur- 
ing combination all spectra were rebinned to a 2.5kms~' pixel 
scale. Both instruments have FWHM (Full-Width Half-Maximum) 
~7 km s~' spectral resolution, thereby giving ~3 pixel sampling for 
unresolved features. 



2.4 Analysis 

The analysis of the final spectra involves two main approaches. 
Firstly, abundance profile analysis based on Voigt profile fits to the 
various ionic transition^] in the absorber and secondly a compari- 
son of the velocity structure of these absorbers to data and models 
in the literature. The first approach allows us to make a quantitative 
analysis of chemical and dust-depletion uniformity and, in princi- 
ple, chemical enrichment history via models of relative abundances 
of a and Fe-peak elements. Unfortunately, the wavelength coverage 
and SNR of the spectra are not sufficient for a full decompositio n 
of the star-formation history (e.g. Dessauges-Zavadsk v et al ■12007b . 
Nonetheless non-uniformity in the a enrichment and dust depletion 



2 http://www.ucolick.org/~xavier/IDL/index.html 

3 uves popler is available from 

http://www.astronomy.swin.edu.au/~mmurphy/UVES_popler.html 



Trans ition wavelengths and oscillator strengths were taken from 
iMortoijfiool . 
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of the absorbers will provide clues as to the physical environment 
of the absorbing gas (See Sections 13,21 & 14.2b . Furthermore, the 
electron density of the absorbing cloud can be constrained from 
measured limits from the N (Fei) /N (Fen) and N (Mgi) /N (Mgn) 
ratios. The second approach involves study of the velocity profiles 
and widths of each absorber, allowing qualitative assessment of the 
physical models able to explain the velocity structure of the gas, 
such as simple galaxy discs, haloes, inflows or outflows. 



2.4. 1 Column densities and equivalent widths 

We used vPFr£| to fit multiple velocity component Voigt profiles to 
each unsaturated metal ion transition of interest in each absorber. 
The transitions of interest for our analysis were Fen, Mgi, Sin, 
Zn ii, Crn, Mn n, Ti n and Al m. An initial guess for the fit was con- 
structed from Fe n as multiple Fe n transitions at a variety of wave- 
lengths were usually observed, which allows vpfit to constrain the 
velocity structure more accurately. For the remaining transitions fit- 
ted, the redshifts and Doppler parameters, b, of individual compo- 
nents were tied to Fe n (or Mg i /12852 when Fe n was unavailable), 
whilst column densities for each component were allowed to vary 
freely. At the resolution and SNR of the spectra we saw no evi- 
dence for these transitions requiring independent redshifts and b 
parameters. 

The total column densities were reported directly from vpfit, 
which allows for a more accurate estimate of the formal error be- 
cause it alleviates the problem of blending between neighbouring 
components, which produces degeneracies between column densi- 
ties and b parameters for those components (see vpfit documenta- 
tion). Upper limits were taken at 3<x using the error reported from 
vpfit. 

Both the Mg 1/12026 and Znn ,12026 and Znii/12062 and 
Cr H/12062 lines are usually blended due to the velocity width of 
the absorption system. Despite this blending the fits to these lines 
are not degenerate; vpfit simultaneously fits these lines, along with 
unblended lines such as Crii/12056 and Mg 1/12852, using the pre- 
determined velocity structure. The spectra are always of sufficient 
wavelength coverage with high enough SNR for a non-degenerate 
fit to be reached. For example the Znn velocity structure may 
appear uncertain if the absorption profile is sufficiently wide to 
blend together the Zn n /12026 and Mg i A2Q26 profile^ and also 
the Zn ii /12062 and Cr n /12062 profiles. However, since the Cr n ve- 
locity structure is independently determined by Crii/12056 and/or 
Cr H/12066, and since fitting the Mg 1/12852 profile determines the 
Mg i /12026 profile, the Zn n velocity structure can be recovered re- 
liably. 

When all transitions of a given ionic species are saturated (e.g. 
Mg n) the column density of the species is unconstrained by Voigt 
profile fitting. Instead, for Mg n /12796, 2803 we measure the equiv- 
alent width of the profile, which gives an indication of the velocity 
spread in the absorber. The equivalent widths were calculated by 
summing across each pixel, i, in the profile 

W = Y j (l-f,)AA, (1) 

i 

where /; is the normalised flux and A/1, is the width of the ith pixel 



5 vpfit is available from http://www.ast.cam.ac.uk/~rfc/vpfit.html 

6 Note that whilst Crii/12026 also contributes to the flux decrement it is 
generally too weak to confuse the deblending process. 



in rest-frame wavelength space. The associated errors are the for- 
mal errors in the sum and thus do not account for errors associated 
with continuum fitting. 

2.4.2 Measuring velocity profile uniformity 

The apparent optica l depth (AOD) method, as first described by 
ISavage & Sembacbl jl99lh . offered an alternative way of measur- 
ing the column densities of velocity components in an absorber to 
more tradi t ional Voigt profile analyses. The reader is referred to 
Prochaska ( 2003) for a discussion of the merits of each technique. 
Prochaski] 1 20031) conducted a velocity bin based analysis of the 
chemical uniformity within DLAs by calculating the column den- 
sity in each velocity bin using the AOD method. Here we will use 
the same technique to measure variation across the velocity profile 
of the absorber in both dust depletion (via Zn and Cr) and chemical 
uniformity (via Fe and Si) where possible. In particular we hope 
to show whether the dust depletion profiles are relatively uniform, 
or whether the dust is concentrated in one or two velocity compo- 
nents as this ma y indicate the presen ce of molecular hydrogen in 
the absorber (e.g.lLedoux et alj2003h . 

|ProchaskaU2003h use unsaturated, unblended transitions such 
as Fe and Si and convert the flux pixels to optical depth pixels, bin- 
ning the result in velocity space. These bins are not components 
in the traditional sense, but just velocity bins across the profile, 
as such they may contain one or more real velocity components. 
Nonetheless, if there is any non-uniformity across the absorption 
pro file, it will still be reflected in the constructed velocity bins. Un- 
like |Prochask3 f2003) we will be applying this technique to blended 
transitions such as Znn and Cm. We circumvent the problem of 
line blending by binning the Voigt profile fits to the data, rather 
than the spectra themselves. The error in each pixel is still drawn 
from the actual spectral error array. This technique utilises the fact 
that we have multiple transitions of each blended element, which 
a straight-forward analysis of the data cannot account for. Using 
the error array corresponding to the flux in each pixel, rather than 
adapting the vpfit code to return the error associated with each pixel 
in the fit, will not underestimate the real error too greatly. This as- 
sertion can be verified by applying the technique to transitions such 
as Fe ii and Si n which do not suffer from blending. A further ad- 
vantage of using the Voigt profile fits to the spectra is that taking 
account of multiple transitions at different wavelengths makes the 
subsequent analysis more resilient to changes in bin size. That is, 
the effect of Root Mean Square (RMS) fluctuations in the data of 
an individual transition is reduced. 

We will use 20km s bins for our analysis, corresponding to 
eight 2.5 km s" 1 pixels per velocity bin; this gives reasonable SNR 
in each bin, whilst still providing enough bins to discuss the ab- 
sorbers' depletion properties. In practice the size of the bin does 
not alter our conclusions as the component-to-component dust de- 
pletions are resilient to changes in bin size; binning mainly effects 
the error bars in our analysis. We also performed the same analysis 
with the velocity bins offset by half a bin from the initial analysis 
to test whether the results were resilient to shifts in velocity space. 
We found this made no difference to our conclusions. 



2.4.3 Electron densities 

The electron density of the gas can be constrained by the relative 
column densities of the neutral and s ingly ionised transitions of 
an element (e.g. lProchaska et al]|2006i) . For instance, under steady 



Ca ii absorber environments 5 



state conditions and photoionization equilibrium the balance of 
photoionisation and recombination for Mg i and Mg n gives 

= AT (Mgi) T(Mgi) 
e N (Mgn) a (Mgn, T) 1 ' 

where N are the measured column densities of the named tran- 
sitions, a (Mgn, T) is the total recombination coefficient (radia- 
tive and dielectric) dependent on temperature, T, and T(Mgi) = 
Const, x tr ph (Mgi) x G/Gq represents the photoionisation rate of 
Mg i. (T p h (Mgi) is cross-section to ionisation integrated over the in- 
cident radiation fi eld and Go = 2.72 x 10~ 3 erg cm' 2 s" 1 is 1.7x 
Habing's constant dHabindl 19681 [Gondhalekar et al.lll980l) . where 
the factor of 1.7 is required to normalise G to the value of the 
Milky Way UV radiation field strength. A similar relationship 
can be derived for Fei & Fen. Thus observations of the ratio of 
the column densities of a neutral and singly ionised species can 
constrain the electron density, n e , for a give n strength of radia- 
tion field, G/G . See IProchaska et alj J200^) Section 4.2 & Ap- 
pendix for details of the calculations. For our analysis we use the 
Fe 1,12484, Fe n 442260, 2374, Mg 1/12852 and Mg n ,12803 transi- 
tions. Note that for the ionisation of Fei, Hi charge transfer can 
do minate the transition, depe ndent on the exact ratio of «hi to 
n e dKingdon & Ferlan3 1 19991) . If Hi charge transfer does domi- 
nate then equation [2] will underestimate the true electron density. 
The results and conclusions drawn from observed limits on elec- 
tron densities are dominated by measurements of Mg i and Mg n 
(See Fig.|5]l- Thus, the increased uncertainty in measurements of n e 
based on Fe due to Hi charge transfer will not effect any conclu- 
sions. 

2.4.4 Constraints on absorber environment from velocity profiles 

There have been extensive studies of the velocity profile structure in 
both DLAs and Mg n absorbers. Early work focused on modeling 
DLAs as rotating discs of gas extending out from the associated 
galaxy dSteidel et al] |2002). Later, more complex models includ- 
ing ho t and cold rotational components, haloes and spherica l accre- 
tion jProchaska & Wolfell997l : lcharlton & Churchilll 19981) or out- 
flows, such as su per-bubbles expa nding out into the Inter-Galactic 
Medium (IGM) (Bon d et all 1200 ll) were called upon to explain 
the increasingly complex Mg n absorption profiles being observed 
and their sometimes highly symmetric velocity profiles. Given the 
paucity of strong Ca n absorbers we do not have a statistical sample 
of absorbers required for a quantitative analysis, therefore we will 
restrict ourselves to a qualitative comparison with models from the 
literature. It will, nonetheless, be possible to place constraints on 
the physical environment of these absorbers such as whether they 
can be reproduced using simple disc models, outflows, etc. 



3 RESULTS 

Voigt profile fits for each absorber are shown in Appendix[A] whilst 
the unfitted Mg n and Ca n profiles are presented in Appendix [B] 

The measured column densities for all the absorbers in the 
sample are given in Table [3] whilst measured equivalent widths 
of Mgn and Can and the velocity width, as measured by AV90 
dProchaska & W olfe 1997), are given in Table [4] The Can transi- 
tion in J1203+1028 is contaminated by sky absorption, but because 
the absorption features are not flat bottomed (i.e. not all the flux is 
absorbed) we are able to make a correction to the W^ 934 measure- 
ment for the sky absorption as follows. The first step was to find 



Table 4. Equivalent width and Avgo measurements from high resolution 
spectra for each absorber in this work. The lcr error in Av'90 is < 10 km s , 
which is dominated by uncertainty in continuum fitting. Where no high res- 
olution data exists the measurement from SDSS spectra are given in paren- 
theses when available. 





Rest-frame equivalent width [A] 




Object 


Mgn 42796 


Can 43934 


AV90 [km s l ] 


J0004-4157 


4.424 ± 0.003 


0.966 ± 0.005 


400 


J0256+0110 


3.27 ± 0.01 


0.39 ± 0.04 1 


350 


J0334-0711 


3.49-4.65 2 


(0.31 ± 0.05) 


410 


J0407-4410 


1.668 ± 0.002 


0.282 ± 0.009 


200 


J05 17-44 10 


2.336 ± 0.002 


0.346 ± 0.002 


470 


J0830+2410 




0.284 ± 0.007 


190 


J0846+0529 


2.22 ± 0.03 


0.50 ± 0.03 


160 


J0953+0801 


0.468 ± 0.0O8 3 


0.34 ± 0.02 


50 


J1005+1157 


2.58 ± 0.02 


(0.78 ± 0.08) 


210 


J1028-0100 


(1.47 ±0.08) 


0.29 ± 0.02 


100 


Jl 107+0048 


2.766 ± 0.009 


0.34 ±0.01' 


190 


Jl 129+0204 


2.11 ±0.01 


0.69 ± 0.03 


170 


Jl 130-1449 


1.794 ±0.007 


0.374 ± 0.005 


120 


J1203+1028 


2.58 ± 0.02 


0.630 ± 0.006 4 


270 


J1211+1030 




0.24 ± 0.01 


330 


J1232-0224 


2.071 ±0.005 


0.203 ± 0.007 


120 


J1323-0021 


2.156 ±0.009 


0.82 ± 0.02 1 


120 


J1430+0149 


2.898 ± 0.009 


0.30 ±0.01 


220 


J2328+0022 


(1.98 ±0.07) 


0.250 ± 0.008 


200 



1 Reported measurement from Nestor et al. (2008) 

2 Velocity structure too complex to separate Mg 11 42796 and Mg 11 42803. 
Therefore we measured the total equivalent width over both lines and as- 
sumed doublet ratios of 1 : 1 and 2: 1 to calculate the range of possible Wg 196 . 

3 W 2796 = 0.88 ± 0.05 A when measured in the SDSS DR4 spectrum of this 
quasar, thus it is included in the sample. 

4 Corrected for sky absorption contamination, as discussed in SectionfJ] the 
associated error only includes the formal statistical error and is therefore an 
underestimate. 

exposures with uncontaminated sky absorption of similar air-mass, 
date and time to the J1203+1028 exposures. For each uncontam- 
inated sky exposure the equivalent width of absorption due to the 
sky, VK sky , was measured. We then apply a correction to the match- 
ing exposure from J1203+1028 such that, W 3934 = W 3934 - W Ay , 
where W 3 ^ 34 is measured across the same wavelength range as W sky . 
The errors are propagated formally, so will not include any system- 
atic error such as that due to continuum fitting. 

3.1 Depletions 

3.1.1 Total depletions 

The total column densities derived for Znn and Cm in Section [3] 
were used to calculate the depletion level in each of the absorbers. 
In high Nm systems self-shielding ensures that the metallic ele- 
ments considered here are predominantly in the lowest ionization 
state with an ionization potential less than 13.6 eV. Thus we may 
assume that the ratio of A'cni/Wznii ~ Na/Nzn, i- e - no correction 
for ionization is required. Measured [Cr/Zn] ratios are given in Ta- 
ble [5] Fig. [T] shows the depletion in each of the systems versus 
W^ 934 for each absorber. Included in the p lot are the DLAs with 
known Zn n , Cm and Ca 11 absorption from Nestor et al. I d2008l) and 
iKhare et al.1 d2004l) . For comparison we sho w the me an values of 
[Cr/Zn] vs. W 3934 measured statistically by |WHP06| from SDSS 
quasar spectra with Can absorbers. The high W^ 934 regime re- 
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Table 3. The ionic total column densities measured for each absorber in this work. Upper limits are quoted at 3<x significance. Systematic errors due to 
continuum level and component fitting are estimated at 0.03 dex. Quoted errors include random and systematic contributions added in quadrature. 



logiV [atom cm 2 ] 



Object 


Can 


Fen 


Mgi 


Znn 


Cm 


Sin 


Aim 


Mmi 


Tin 


J0334 




14.94 ± 0.04 


12.84 ±0.07 


12.58 ± 0.09 


13.19 ±0.08 






12.65 ± 0.04 


11.9 ±0.1 


J0517 


12.74 ± 0.03 


14.31 ±0.03' 


12.74 ± 0.03 


12.22 ± 0.04 


12.53 ±0.06 


14.77 ± 0.03 2 


12.76 ± 0.03 


12.21 ± 0.03 


<10.6 


J0846 


13.06 ± 0.06 


15.21 ± 0.04 


12.96 ± 0.05 


<12.8 


<13.4 






13.21 ± 0.06 


13.00 ± 0.05 


J0953 


13.57 ±0.04 


15.09 ± 0.04 


13.24 ±0.03 


12.25 ± 0.06 


13.17 ±0.04 


15.61 ± 0.07 


15.31 ± 0.03 


13.0 ±0.7 


12.1 ±0.1 


J 1005 




15.30 ± 0.05 


13.10 ±0.03 


<12.8 


<13.3 




<13.5 


13.03 ±0.05 


<12.4 


J1107 




15.52 ±0.03 


13.00 ± 0.07 


13.23 ±0.04 


13.83 ± 0.09 






13.37 ±0.04 


13.12 ±0.03 


J1129 


13.11 ±0.04 


15.35 ± 0.03 


13.14 ±0.09 3 


12.80 ± 0.05 


13.66 ± 0.04 


16.0 ±0.5 


13.81 ±0.09 


13.13 ±0.04 


12.79 ± 0.04 


J 1203 


13.05 ± 0.03 


14.83 ± 0.03 


12.99 ± 0.03 


12.63 ± 0.06 


13.38 ± 0.05 






12.77 ± 0.05 


<11.2 


J1430 


12.85 ±0.04 




13.07 ± 0.05 


12.94 ± 0.04 


13.45 ± 0.04 


15.75 ± 0.03 


13.72 ± 0.03 




12.79 ± 0.04 



' Fen measured over full ~ 800kms~' range of absorber. A^eii = 14.25 ± 0.03 over ~ 100kms~' region which other elements could be measured. 

2 Sin constrained by simultaneous fit to Siii/11526, 1808 and intervening Lya forest. 

3 Mg i /12852 not covered by spectrum. Mg i column density derived from residual from Zn n and Cr n profile fits. 
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Figure 1. [Cr/Zn] vs. W 3 for each of the absorbers. The black diamonds 
show [Cr/Zn] measured in the UVES or HIRES spectra vs. W 3934 measured 
from the SDSS spectra. An indication of the typical error is shown at the 
centre-to p of the p lot. The dashed grey regions show the average relation 
found by WHP06. Other measurements of [Cr/Zn] fo r Can absorbers in 
the literature are shown as black squares and asterisks INestor et al . 2008; 
Khare et alj2004l) . The horizontal dashed black line shows the level of solar 
depletion. The vertical dashed grey line shows the W 3 ,' 34 equivalent width 
limit of our sample selection. 



Figure 2. [Cr/Zn] vs. z abs for the DLAs from lAkerman et alj fe005r) and 
each of the absorbers. The black diamonds show [Cr/Zn] measured in the 
UVES or HIRES spectra. The blue crosses show the DLA measurements 
of [Cr/Zn] whilst the blue triangles show upper limits to the depletion. The 
dashed black line shows the level of solar depletion. 

warm halo phase in either the MW, LMC or SMC, the [Si/Zn] ratios 
are ~ 0.5 dex higher than typically observed locally. 



mains relatively unexplored for measurements from high resolution 
spectroscopy. Similarly, Fig. [2] shows [Cr/Zn] vs. redshift for the 
Can absorbers wi th the compilat i on of measurements of [Cr/Zn] 
in DLAs given by lAkerman et al.l J2005h included for comparison. 
We note that the measurements from the observed Can absorber 
sample seem consistent with those observed in other DLAs, being 
neither particularly dusty nor dust-free. Performing a Kolmogorov- 
Smirnov (KS) test of the measured [Cr/Zn] ratios indicates an 87 
per cent likelihood that the Ca n and DLA samples are drawn from 
the same underlying population. 

Table [5] also shows several abundance ratios sensitive to dust 
depletion. It is evident that the depletion properties of these ab- 
sorbers diff er slightly from loca l measurements of the MW, LMC 
and SMC ISavage & Sembachl 1 19961 : IWeltv et alj 1 19991 1200 ll) : 
whilst the [Zn/Fe] ratios in these Ca n absorbers are similar to the 



3.1.2 Depletion uniformity 

Using the apparent optical depth method detailed in Section l2.4.2l 
we calculate the variation in [Cr/Zn] across the velocity profile in 
each absorber for which both Znn and Cm are detectable. The re- 
sults are shown in Fig. [3] 

There is little evidence for any significant (> 0.3 dex) devia- 
tion from uniform depletion across the velocity profile in each of 
the absorbers, even for those absorbers with reasonably broad ve- 
locity profiles (> 50kms~'; the exception is J05 17-4410, to be 
discussed later). It is not necessarily surprising that an absorber 
with complex velocity structure might appear well mixed because 
velocity components do necessarily not reflect the spatial structure 
of the absorbing gas. 0.3 dex variation has 3<r significance at the 
SNRs of the spectra, thus the results in Fig.[3]restricts the range of 
dust depletion models which may apply to these Can absorbers. In 
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Table 5. A table of abundance ratios in the absorber sample, relative to solar, traditionally used for assessing 
dust-depletion and chemical enrichment history. Solar abundances are taken from Lodders (2003). 



Object 


[Cr/Zn] 


[Mn/Fe] 


[Ti/Fe] 


[Zn/Fe] 


[Si/Fe] 


[Si/Zn] 


[Si/Ti] 


J0334-0711 


-0.4 ±0.1 


-0.33 ± 0.06 


-0.5 ±0.1 


0.5 ±0.1 








J0517-4410 1 


-0.71 ±0.07 


-0.08 ± 0.04 


< -0.70 


0.81 ±0.05 


0.44 ± 0.04 


-0.36 ± 0.05 


> 1.14 


J0846+0529 




-0.02 ± 0.07 


0.35 ± 0.06 


< 0.62 








J0953+0801 


-0.10 ±0.08 


-0.2 ± 0.7 


-0.5 ±0.1 


0.00 ± 0.07 


0.45 ± 0.08 


0.45 ± 0.09 


0.9 ± 0.1 


J1005+1157 




-0.30 ± 0.06 


< -0.39 


< 0.32 








Jl 107+0048 


-0.4 ±0.1 


-0.17 ±0.06 


0.17 ± 0.05 


0.53 ± 0.06 








Jl 129+0204 


-0.16 ±0.06 


-0.26 ± 0.05 


-0.01 ± 0.05 


0.29 ± 0.06 


0.5 ± 0.5 


0.2 ± 0.5 


0.5 ± 0.5 


J1203+1028 


-0.27 ± 0.06 


-0.08 ± 0.06 


< -1.08 


0.64 ± 0.07 








J1430+0149 


-0.51 ±0.06 










-0.11 ±0.05 


0.34 ± 0.05 



Abundance ratios for this absorber based only on regions in velocity space where both transitions are de- 
tectable; in particular the Fe n transitions cover a much broader region of velocity space than most other transi- 
tions. 



many cases it is possible to restrict the variation to < 0.2 dex varia- 
tion in [Cr/Zn], i.e. there is no variation with significance > 2a for 
most absorbers. See Section HTTI for further discussion. We note that 
the velocity width of J0953+0801 is too narrow for any assessment 
of variation across its profile with the bin size necessary for the 
analysis (See Fig. |3(b)[ >. In addition, the Mg 1/12026 contribution to 
the flux decrement in the region of Znii/12026 for J1129+0204 is 
uncertain. The uncertainty arises because the Mg i contribution to 
the flux decrement in the blended part of the spectrum could not 
be constrained from the Mgi/!2852 line, which is not observed. 
Instead the Mg i /12026 line is constrained as the residual from fit- 
ting the Zn ii and Cr n lines alone. As can be seen from Fig. [3(d)] 
the variation in [Cr/Zn] traces regions with the most contamina- 
tion from Mgi. Hence, the possible signature of variation across 
the Jl 129+0204 profile in Fig. |3(d)| cannot be considered robust. 

Fig.[3]presents two different velocity binnings, each with the 
same bin size but offset from each other by 1 km s~ 1 . The similar- 
ity between the two binnings shows that the analysis is resilient to 
shifts in the bin centres. Moderate changes, of order 10-20kms~', 
to the size of bins also has little effect on the results, other than al- 
tering the magnitude of the associated error bars. 

J05 17-4410 has substantial variation in [Cr/Zn] across its 
~ 800kms~' profile. J05 17-4410 is also a known sub-DLA with 
a high UV ionising field JOuast et aljp200^ , thus much of the ob- 
served variation is likely due to differential ionisation of the gas, 
rather than being purely due to differential depletion of the gas 
onto dust grains. Most transitions are only detected in the central 
~ 100 kms~', the strongest part of the absorber, thus we will con- 
centrate our analysis in this region (See Fig, |3(c)fr . The advantage 
of restricting the analysis is that the gas may exhibit less variation 
due to ionisation on smaller scales and over this velocity range the 
absorption profile is relatively uniform, except for a few velocity 
bins near -25kms~', which exhibit ~ 0.3 dex (3cr significant) de- 
viation. 



3.2 Chemical history 

3.2.1 Enrichment history 

We have detected both ff-capture elements and Fe-peak elements 
in high resolution spectra of these Can absorbers, a elements 



such as Si, Ca and T0 are created by type II supernovae on a 
timescale of < 2 x 10 7 yr after a star formation event, whereas 
the Fe-peak elements are created by type la supernovae, which oc- 
cur on a timescale of 10 8 - 10 9 yr. Thus, the ratio of a to Fe-peak 
elements can reveal something about the en richment history and 
previous star-formation in the a bsorber (e.g. IWheeler et alj|l98^ : 
iDessauges-Zavadskv et alj|2007h - Unfortunately the same element 
ratios are also affected by dust-depletion and disentangling the two 
effects is not trivial. However, if one observes enough different 
transitions it is possible to begin to separate the two effects (e.g. 
iHerbert-Fort et al]2006l) . Without a measurement of the metallicity 
of the absorber, a full deconstruction of the star-formation history 
of the absorber is not possib le. Therefore, we follow the example 
of IHerbert-Fort etHll feOOcjj) , presenting various abundance ratios 
which are instrumental in constraining the effects of dust-depletion 
and enrichment history in Table[5] 

The abundance ratios measured in the sample of Can 
absorbers are general ly consistent with the DLA population 
jProchaska et alj|2003l) . Furthermore, up to 0.33 dex variation is 
to be expected purely from differential dust depletion in the sam- 
ple, particularly for elements strongly affected by dust suc h as Ti 
dSavage & Sembachll 19961 : IDessauges-Zavadskv et al .1 12006 ) . 

The Metal - Strong DLAs (MSDLAs) studied by 
IHerbert-Fort et al] fc006l) represent an unusual subset of DLAs and 
so we consider whether Ca n absorbers overlap at all with these rare 
systems, as similarities between the populations may provide clues 
as to the physical origin of th e absorbers. To cla s sify a n absorber 
as a MSDLA, as defined by IHerbert-Fort et all J20061) . requires 
Af(Znii) > 13.15 or N (Sin) ^ 15.95. Under this definition only 
two of the absorbers (11107+0 048 and Jl 129+020 4) marginally 
classify as MSDLAs. Similarly, iNestor et alj d2008h measure five 
of nine Can absorbers to be MSDLAs, although they are also all 
marginal cases. Therefore, with the caveat that the highest Wg 934 
systems remain unprobed, the two populations are obviously not 
identical. Rather, it seems that the tail end of the Ca n metal column 
density distribution overlaps with the definition of MSDLAs. If 
some strong Can absorbers arise, not because of sputtering of 
dust grains or weak ionisation backgrounds (see introductions of 



7 Ti is not strictly speaking an a element, but it exh ibits abundance pat- 
terns similar to othe r a-elements in Galactic stars I Edvardss on et al.1119951 ; 
iFrancois et alj2004h . 
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J0334 -0711, z = 0.599, [Cr/Zn] = -0.412 ± 0.109 



J095J+0801, z = 1.023, [Cr/Zn] = -0.101 ± 0.054 
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(a) The depletion profile across the z, = 0.59760 absorber towards J0334- 
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(b) The depletion profile across the z = 1.02316 absorber towards J0953+0801. 
The velocity profile for this absorber is too narrow to measure variation in 
[Cr/Zn]. 



J0517-4410, z = 1.151, [Cr/Zn] = -0.710 ± 0.050 



J1 129+0204, z = 0.965, [Cr/Zn] = -0.157 ± 0.045 
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(c) The depletion profile across the z = 1.14955 absorber towards J0517-4410. 
There is evidence for significant variation, ~ 0.3 dex (3<x), in th e absorber's 
[Cr/Zn] profile, consistent with findings of previous works (e.g. Ipuast et al] 
l2008h . 
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(d) The depletion profile across the z = 0.96497 absorber towards Jl 129+0204. 
Mg i /12852 was not observed so the profile fits for Zn n and Cr n are less certain, 
making it difficult to assess profile uniformity. 



Figure 3. Each top panel shows the Crii/12056 transition, each middle panel shows the blended Znii/12026 and Mg 1/12026 transitions. The flux is shown as 
black histograms, the continuum as a dashed line (blue) and the fit to the Cm or Znn transition is shown as a solid line (blue). The dashed line (red) through 
the data shows the combined fit to Mg i and Zn n. Each bottom panel shows the value of [Cr/Zn] across the profile. The short-dashed line (blue) represents the 
solar depletion level whilst the dot-dashed line (red) shows the total depletion measured in this absorber as referenced in Tableland Section [3.1.1l The black 
points are one set of 20kms~' bins and the blue crosses are a different set of similarly sized bins, offset by 10kms~' from the first. The value of^ 2 given in 
the top left-hand corner of the panel is the^ 2 per degree of freedom for the fit of the red dot-dashed line to the black points. That is, x 2 represents how well 
the profile is fit by a single depletion value. 



WHPOI IZvch et all 120071) but simply because of large column 
densities of metals, this would be expected. 



3.2.2 Chemical uniformity 

Given the relative uniformity of the depletion profiles, as presented 
in Section [3.1.2l any variation in the [Fe/a] ratio (where a is any 
a-capture element, such as Si, Ti or Ca) should be indicative of 
chemical non-uniformity in the absorbers, rather than variation in 
dust depletion. For the absorbers in our sample where we can mea- 
sure [Fe/o-] the results are presented in Fig. [4] Even without metal- 
licity measures it is possible to see abundance variations across the 



profile, which may indicate that different components of the gas 
have different enrichment histories. In the case that there is no sim- 
ilar variation in the [Cr/Zn] ratio at the same velocities then the 
variation is most likely due to or-enhancement. 

The use of the unsaturated, unblended Fen and Sin lines in 
this analysis also allows us to test the validity of using the fit- 
ted absorption profile, rather than the data themselves, to calcu- 
late the blended [Cr/Zn] depletion profiles. We observe no signifi- 
cant « 0.1 dex difference to the derived [Fe/Si] ratios either when 
analysing the Voigt profiles fits or when analysing the data directly. 

Only three absorbers have sufficient wavelength coverage 
to study their [Fe/Si] ratios and we comment on each in turn. 
J0953+0801 is not broad enough in velocity space to show signif- 
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J1 107+0048, z = 0.740, [Cr/Zn] = -0.350 ± 0.145 
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(e) The depletion profile across the z = 0.74030 absorber towards Jl 107+0048 
J1203+102B, z = 0.746, [Cr/Zn] = -0.274 ± 0.064 
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(f) The depletion profile across the z = 0.74630 absorber towards J1203+1028. There is a blend with a z = 1.322 C iv absorber at -1 10 and -130kms~'. 



J 1430+0149, z = 1.242, [Cr/Zn] = -0.512 ± 0.039 
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(g) The depletion profile across the z = 1.24180 absorber towards J1430+0149. There is a blend with a z = 1.933 C iv absorber at — 120 km s - '. 

Figure 3 - continued 



icant variation across its profile. Jl 129+0204 shows a ~ 0.4 dex 
(3cr) deficit of Si toward the left hand side of its profile (See 
Fig. [4(b)}, however it has questionable significance; see discussion 
in Section|4~2l J0517-4410 also shows > 0.8 dex (8cr) variation in 
its [Fe/Si] profile, but this variation is mirrored by its [Cr/Zn] pro- 



file (See Fig. |3(c)} , which makes it likely the variation is a further 
signature of differential dust depletion, rather than a enhancement. 
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J0517-4410, l m = 1.151, [Fe/SI] = -0.444 ± 0.008 
1.2 c ■ ■ ■ — 
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(c) The [Fe/Si] profile across the z = 1.14955 absorbers towards J0517-4410. 
Note that the leftmost velocity feature seems to have a different chemical history 
to the rest of the profile. Given the [Cr/Zn] profile (see Fig. |3(c)) shows similar 
variation, the [Fe/Si] variation is likely an effect of dust depletion, rather than 
differential enrichment history. 



Figure 4. Each top panel shows the Feii/i2250 transition, each middle panel shows the Siii/11808 transition. The flux is shown as black histograms, the 
continuum as a dashed line (blue) and the fit to the Fe n or Si n transition is shown as a solid line (blue). Each bottom panel shows the value of [Fe/Si] across 
the profile. The short-dashed line (blue) represents the solar abundance ratio whilst the dot-dashed line (red) shows the total abundance ratio measured in this 
absorber as referenced in Tableland Section POl The black points are one set of 20 km s~ 1 bins and the blue crosses are a different set of similarly sized bins, 
offset by 10 km s from the first. The value of \ 2 given in the top left-hand corner of the panel is the x 1 per degree of freedom for the fit of the red dot-dashed 
line to the black points. That is, x 1 represents how well fit the profile is by a single abundance ratio value. This absorption profile is too narrow to assess any 
signature of chemical non-uniformity. 



3.3 Electron densities 

Given the wavelength coverage of the spectra in our sample, 
there are two sets of transitions which can be used to con- 
strain the electron density in each of our absorbers: Fei/i2484 & 
Fe ii A/12260, 2374 and Mg i ,12852 & Mg n A2803. It is not possible 
to measure an exact ratio for either set of transitions, just limits. 
Fei is not detected at the SNR of the spectra, thus the measured 
N (Fei)//V (Fen) ratio is an upper limit. Mg n is saturated, thus the 
measured N (Mgi)//V (Mgii) ratio is also an upper limit. Both ratios 
are derived using column densities measured via the AOD method 
(See Table[6]for results). 

Using the AOD method ensures that the column densities for 
each transition are measured over the same velocity range, provid- 



ing robust limits. For saturated pixels a minimum normalised flux is 
assumed of 0.05 or cr(flux) /5, whichever is greater; this provides 
a conservative lower limit to the optical depth. For each absorber 
we derive the product n e (G/Go) 1 using a conservative estimate 
for the temperature, T = 8000 K, which maximises the recombina- 
tion rates. Fig.|5]shows n e (G/Go)~' versus W^ 934 for each absorber. 
We take the limit from Fe or Mg, whichever provides the greatest 
constraint on n e in each case. Fig. [5] indicates that n e < 0.1 ctrT 3 
in Can absorbers if we assume G/Gq ~ 1, as in the Milky Way. 
Thus, Can absorbers are not comprised of extremely dense mate- 
rial since iiu ~ latomctrT 3 , unless the gas is extremely neutral 
(n H +/"H ^ 0.01). Here we assume that the free electrons balance 
the ionized hydrogen, n e / (rc H + / n n) = "h- 
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Figure 6. The velocity structure in each of the absorbers, as represented by an unsaturated, unblended absorption line, identified in the top right-hand corner 
of each plot. The absorber's quasar is identified in the top left-hand corner of each plot. The black histograms are the data and the red histograms are the error, 
whilst the continuum is the green short-dashed line. For the highly complex absorbers towards J0334-071 1 and J05 17-4410 we show a second transition with 
a different /-value in grey to reveal more of the structure. 



3.4 Absorber environment 



As detailed in Section [2.4.4l modelling absorbers as rotation discs 
or outflows, etc. allows us to predict the velocity structure of such 
gas. Thus, it is possible to constrain the physical environment of 
the absorber by ruling out mechanisms which cannot reproduce 
observed velocity structures. The velocity profiles of the Can ab- 
sorbers are shown in Fig. [6] Here we present any unsaturated tran- 
sition with high SNR, such as Mg i ,42852 or Fe n /12374 from each 
absorber for comparison. There are a wide range of different pro- 
files associated with Can absorbers, so it is unlikely that a single 



physical process gives rise to the strong absorption; see Section l4~4l 
for further discussion. 



4 DISCUSSION 
4.1 Depletions 

The Can absorbers observed appear to have very similar deple- 
tion properties to the general DLA population (See Fig. [3. Whilst 
this similarity may at first seem contrary to the results of WHP06, 



12 B. J. Zych et al. 



Table 6. The column density ratios of neutral-to-singly ionised species in 
the Can absorber sample. 
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Figure 5. Upper limits on the product of the electron density n e and the 
near-UV radiation field G scaled to 1.7x Habing's constant Go = 2.72 X 
10~ 3 er g cm -2 s -1 , which roughly corresponds to the Galactic far-UV in- 
tensity (Habing 1968; Gondhalekar et al. 1980). The limits are established 
by our observed upper limits to the N (Fei)/N (Fen) and N (Mgi)/N (Mgn) 
ratios under the conservative assumption that the electron temperature is 
T e = 8000K. Limits based on Mg are black, whilst those based on Fe are 
grey. The results indicate the Ca n absorbers are not comprised of extremely 
dense material. 



who find Can absorbers to be amongst the most dusty absorbers 
known, Fig. Q] shows the results are cons istent. Ou r observations 
only probed the lower W^' 34 regime, which lWHPOrj find have much 
smaller levels of depletion. The distribution of po ints in Fi g.QJmay 
suggest that the average depletions measured by|wHP06 define a 
lower envelope to the dust depletions of Can absorbers, rather than 
a simple linear correlation. That is, low W^ 934 absorbers may com- 
prise both weakly and strongly dust depleted systems, while high 
Wq 934 absorbers may all be strongly dust depleted. Such a situation 
may point to two distinct mechanisms producing Can absorbers. 
The seven Can absorbers for which we were able to measure 



the dust depletion properties across the velocity profile have varia- 
tion < 0.3 dex (corresponding to 3tr) in [Cr/Zn]. Once again, these 
dust depletion properti es are consistent with results from the gen- 
eral DLA population (Prochaska 2003). Given the uniformity of 
the dust depletion it is unlikely that any of these absorbers harbour 
molecular hydrogen since we would expec t the dust depletion to be 
concentrated into one or two components dLedoux et alj2003h . 

J0517-4410 is exceptional in its depletion variation, evident 
in both the [Cr/Zn] and [Fe/Si] profiles (See Figs. |3(c)1 & gg). 
Molecular hydrogen h as been detected in this absorber at z = 
1.15079 and 1.15085 dReimers et ail [2003), which matches with 
the velocity bins exhibiting increased dust depletion. Note that nar- 
row C i absorption lines have also been detected at the same red- 
shifts dLevshakov et"ai]|2006l) . It is possible that H2 traces these 
neutral components, rather than dust as measured using Znn and 
Crn. Given that the AOD velocity bins would smear out the de- 
pletion signature, particularly in very narrow, cold components, the 
accuracy with which [Cr/Zn] traces H 2 or Ci cannot be assessed 
more precisely here. The abundance ratios from our analysis of 
J05 17 -4410 agree with the more detailed analysis of lOuast et al.l 
d2008h . The abundance ratios relative to Zn in this absorber match 
very closely the SMC abundances towards Sk 155 (gro up A) which 
are co nsistent with MW halo cloud abundance ratios ( iWeltv et"aT1 
l200ll) . except for [Si/Zn] ~ -0.36 dex and [Ti/Zn] < -1.9, which 
are somewhat depleted. The significance of these differences is dif- 
ficult to determine because analysis of the [Si/S] ratio along the line 
of sight to Sk 155 has shown that Si i s more depleted toward Sk 155 
than was apparent dSofia et alJibOo!) . 

The lack of evidence for molecular hydrogen from non- 
uniform d ust depletion is seem ingly at odds with the recent re- 
sults from iNestor et alj d2008l) who use photo-ionisation model- 
ing to show that Can absorbers probably have total gas densi- 
ties > 1 atom cnr 3 - the re gime where we woul d expect molec- 
ular hydrogen to be present dSrianand et al.ll2005l) . There is tenta- 
tive evidence that stronger W„ 934 systems should have larger total 
gas densities dNestor et alJ|2008T) , therefore the results presented 
here may just reflect the fact we have only probed the lower W^' 34 
regime and that it is the stronger WZ 93 * systems which are associ- 
ated with molecular hydrogen. Another possibility is that we need 
to rea ssess some of the as sumptions in the photo-ionisation model- 
ing of lNestor et alJd2008l) . For instance Ne stor et alj d2008f) assume 
that the UV radiation field in the vicinity of their absorbers is simi- 
lar to that of the solar neighbourhood; perhaps a more general range 
of UV fields needs to be considered (see Section |4~3l . 



4.2 Chemical history 

Detailed analysis of the star formation history in Can absorbers is 
not appropriate given the limitations of the spectra; a wider range 
of transitions, both Fe-peak and a elements, spanning a range of 
ionisation potentials, as well as observations of the H 1 column den- 
sity to measure metallicities are required to conduct a full anal- 
ysis of the star formation history. It is, nonetheless, possible to 
use the [a/Fe] ratio across the velocity profile to infer whether 
the gas in the absorber has multiple star-formation histories. Only 
two absorbers, Jl 129+0204 and J05 17-4410, showed any varia- 
tion in [Si/Fe] > 0.3 dex (3cr significance) across their profiles (see 
Figs. |4(b)| & |4(c)| l. Given that dust depletion can also cause vari- 
ation in [Si/Fe] as well as differential enrichment history it is im- 
portant to compare the [Si/Fe] and [Cr/Zn] profiles to break this 
degeneracy. 

Recall that the dust depletion profile of Jl 129+0204 
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(Fig. |3(d)} has ~ lcr significant variations at the velocities where 
the [Fe/Si] profile has ~ 3tr significant deviation at -25kms~'. 
It appears that the components from -100 to -25kms~' may be 
offset in [Fe/Si] from the components from -5 to 110kms~'. The 
apparent deviation is not significant given the a posteriori nature 
of the statistic. Only higher SNR data where all the individual bins 
between -100 and -25kms~' are > 3<x from the mean can con- 
firm whether this deviation is real. Such a deviation would be of 
interest because it is indicative of the gas having different star- 
formation histories. Possible explanations would include that the 
absorbing gas arises from the merging of two galaxies with differ- 
ent star-formation histories. Or we could be observing chemical in- 
homogeneities within a single galaxy; certainly our own Milky Way 
is obs erved to have variations in a-enhancement (e.g. lCunha et al.l 
120071) . Alternatively, we could be witnessing the building blocks of 
a galaxy coming together at z ~ 1. Such assessment must wait for 
higher SNR data however. 

The [Fe/Si] profile for J0517-4410 has > 0.8 dex (8tr) devia- 
tion near -25 km s from an otherwise relatively uniform profile. 
Note, however, that a similar ~ 0.3 dex (3<x) deviation is seen in 
the [Cr/Zn] profile. Thus it is likely that this is purely a signature 
of dust depletion, rather than cr-enhancement in the absorber. Note 
that J05 17-4410 i s a sub-DLA with a strong ionising background 
dOuast et al.ll2008h thus it is possible that AWi/A'siii £ AW_/Vsi and 
some of the observed deviation in [Fe/Si] and [Cr/Zn] could be due 
to variable ionisation. In addition, the trend in [Fe/Si] with N (Sin) 
follows what one would expect if photo-ionisation were causing the 
variation. Nonetheless, across the region we have analysed here, the 
fraction of highl y ionised t o lowly ionised gas appears relatively 
constant dOuast et al.ll2008l on-line material), thus it is likely the 
variation is due to dust depletion, rather than ionisation. This argu- 
ment is strengthened by the fact that molecular hydrogen and Ci 
are detected at the same velocities (See Section |4~0 . 

4.3 Electron densities 

In Section [3~3l we show that the Can absorbers in our sample have 
n e < 0. 1 cm 4 . This relies on the assumption that the radiation field 
in these absorbers is similar to or weaker than the radiation field 
of the Milky Way. Whilst in a few cases the radiation field may be 
greater than this, it seems unlikely to be the case in most absorbers 
as a stronger radiation field would lead to more highly ionized gas, 
thus this remains a robust upper limit. Note that the extragalactic 
UV backgr ound is about lOx weaker in the far UV than the Galac- 
tic value dSuiatha et alj|2008h thus absorbers experiencing strong 
UV radiation fields must be associated with galaxies emitting a lot 
of energy in the far UV. 

The electron densities and gas densi ties derived i n Sec - 
tion 13.31 are within t he range predicted by iNestor et al.l d2008h . 
Srianand et al. ( 2005) find molecular hydrogen to be associated 
with absorbers whose gas densities are n H > 1 atomcirr 3 , which 
is similar to the probable gas densities in the Ca n absorber sam- 
ple derived in Section 1331 «h ~ 1 atomcirr 3 . Fig. [5] provides lit- 
tle evidence for gre ater gas densities at stronger W^ 934 , as found by 
Nestor et al. (2008). However, we can only place limits on n e , rather 
than direct measurements so a detailed comparison is not possi- 
ble. For instance, if stronger Wq 9M absorbers preferentially arise in 
stronger radiation fields for a given ionisation fraction, then Fig. [5] 
would imply greater «h at higher W^ 934 . 

Three absorbers in our sample, those towards J0256+0110, 
Jl 107+0 048 and J1323-002 1, overlap with the absorbers mod- 
elled by INestor etaL | d2008l) . The cloudy models predict spe- 



cific electron densities for each absorber: 0.050, 0.035, and 
0.058 cm- 3 for J0256+0110, Jl 107+0048 and J1323-0021, re- 
spectively (D. B. Nestor, private communicati on). Combining the 
electron densities measured by Nest or et"al] with the limits on 
n e (G/Go) 1 derived in this paper we place lower limits on the radi- 
ation field strengths of 2.6x, 1.6x and 0.5x (G/Go) in the absorbers 
towards J0256+0110, Jl 107+0048 and J1323-0021, respectively. 
Further measurements would be required to comment on whether 
there is a trend with W^ 934 and radiation field strength, though it is 
interesting to note that the weakest field is associated with the ab- 
sorber with the highest W^ 934 of the three. Note that these are only 
rough es timates of the field strength as the electron densities mea- 
sured by INestor et al.l depend on the strength of the radiation field 
used in the cloudy modelling. That is, running the cloudy models 
again with a different value for G/Go will result in a different value 
of n e returned from the model, which in turn will lead to a new 
determination of G/Go- 

Generally Mg places stronger constraints on n e than Fe in 
the spectra and, in any case, all the constraints are li mits, so it is 
difficu lt to compare them with the measurements of INestor et al.l 
(2008). Constraining n e and hence «h more accurately might be 
achieved by spectral observations of the C n* and Si n* fine structure 
lines. However, note that these trans itions constrain n e more tightly 
if the gas is predominantly ionised (Prochaska 1999). Higher ion- 
isation state lines such as Civ, may arise in distinct regions and 
therefore not bear too directly on the Can absorbing gas. 



4.4 Environment 

Velocity structure by itself gives us limited insight into the true 
physical environments hosting Can absorption as we can repro- 
duce the velocity profiles of all absorbers via various models in 
the literature (Recall Section [2.4.4l >. A further consequence of such 
models is that it is possible for multiple physical processes to pro- 
duce similar velocity profiles. It is even possible to take a patho- 
logical line of sight through an outflo w or inflow to reproduce wha t 
looks like a simple disc profile (e.g. Charlton & Churchill 19981) . 
However, it is not possible for a pure disc model to reproduce 
an outflow-like profile because not enough high-velocity com po- 
nents are recovered from the model (e.g. Kacprzak et al. 20081). At 
a cursory glance it seems that the Ca n absorber sample has much 
greater velocity widths than the general DLA population; that is, 
they are less disk-like than general DLAs. Taking a sample of 95 
DLAs with z > 1.6 b ased on observations fr om ESfl HIRES and 
UVES as described in lProchaska et al.l d2003l) we can measure their 
veloci ty widths using Avjo, as described by [Frochaska & Wolfd 
(1997). We find that the median velocity width of such as sam- 
ple is (Av) DLA =i 75kms~', whereas the Can absorber sample pre- 
sented here has a median velocity width, (Av) Can - 200 km s -1 (See 
Table [4] for Av 90 measurements). Although we made no conscious 
decision to select high velocity width absorbers the Can sample 
is, in fact, biased due to the limit of detecting a 0.2 A equivalent 
width feature at 4cr significance in SDSS spectra. The exact ve- 
locity cut this implies depends on the depth of the Can absorption 
and because we do not have this information for the DLA sam- 
ple it is non-trivial to apply a cut to those data. For instance, a 
0.6 A equivalent width feature in a spectrum with SNR~ 20, would 
require Av 90 ~ 120kms _1 to be detectable. Furthermore, if one 
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Figure 7. The meas ured Ai'90 distribution for the DLA sample of 
Proc haska etafl [2003), shown as the black histogram. The AV90 distribu- 
tion for our sample of Ca 11 absorbers is also shown as the hashed-grey his- 
togram. The number of observed Ca 11 absorbers is too small for a meaning- 
ful statistical comparison of the two populations. Given the Wg 934 detection 
limit for Can absorbers implies a cut in AV90, the two absorber populations 
could have the same intrinsic AV90 distributions. The approximate cut im- 
plied in Av'9o at 150 km s is shown as a dotted vertical line (See discussion 
in Section IS) , 



compares the distribution of AV90 measurements of the two sam- 
ples (See Fig. [7} it is evident that the number of DLAs dramatically 
increases at Avqo < 120 km s . Given the Wg 934 detection limit for 
Can absorbers implies a cut in Avgo and the small sample size, the 
two absorber populations could have the same intrinsic Avjo dis- 
tributions. Above 150 kms~', where velocity-selection effects are 
less important, there is a 78 per cent chance the two Av 90 distribu- 
tions are drawn from the same underlying distribution, based on a 
KS-test. 

Note that the selection criterion, Wg 796 ^ 0.5 A implies a ve- 
locity cut on the Can sample. If we conservatively assume that the 
Mg 11 absorption is completely saturated, then the minimum veloc- 
ity width is AV90 ~ 50kms~', for our SDSS selected sample. Our 
Can absorber sample includes systems down to this limit. This se- 
lection effect will not be as important as the selection effect due to 
the Wq 934 limit of any larger sample studied in the future because 
the bulk of DLAs observed have Avjo greater than 50 km s . 

Irrespective of velocity-selection effects, it is clear there are a 
mixture of narrow and wide velocity profiles in the Can absorber 
sample (See Fig[6). The above result is suggestive that strong Can 
absorbers consist of two populations, one which is due to a line-of- 
sight passing at low impact parameter through a quiescent galaxy 
disc and the other due to a line-of-sight passing through a complex 
outflow, merger or galaxy cluster. As a further point, Can absorp- 
tion has been detecte d in many High Velocity Clouds (HVCs ) in 
the Milky Way (e.g. IWeltv etafll 1996t iBen Bekhti et al.ll2008b - If 
one modelled a galaxy as a disc + HVCs, one would still not be 
able to reproduce the high velocity components of the extragalactic 
Can absorbers because the H^' 34 of HVCs are too small. 

Simulations o f DLAs in a cosmological setting (e.g. 
IPontzen et al.ll2008h significantly under- predict the relative num- 
ber of complex absorber profiles compared to simple disc-like pro- 
files, suggesting that we do not fully understand the relative impor- 
tance of different physical mechanisms for producing absorbers. 
Nonetheless, most, if not all, Can absorbers studied here are not 



simple disc-like systems as their Fe 11 absorption line profiles are 
too complex (See Fig. [6}. 



5 CONCLUSIONS 

In this work we present the first high resolution study of Can- 
selected absorption systems, using spectra from the UVES and 
HIRES echelle spectrographs. We divided our analysis into three 
parts, examining the properties of dust, cr-enhancement and envi- 
ronment of the absorbers in turn, finding the following: 

(i) We demonstrate for the first time the feasibility of using 
AOD methods on blended absorption profiles, such as Znii/12062 
and Cr H/12062, by utilising the Voigt profile fits, rather than the 
flux, which enables the assessme nt of velocity p rofile uniformity. 
We find, as previous authors (e.g. IProch aska 2003), that the AOD 
method is robust against binning size and positioning. Furthermore, 
we show there is little difference in applying the AOD to the fluxes 
compared to an absorption profile fit to the spectra, thus allowing us 
to use the method on blended lines. The Voigt profile AOD method 
allows us to combine the power of the more traditional Voigt profile 
fitting and the AOD methods when studying the absorber. 

(ii) From a study of dust depletion across the absorber profile us- 
ing the Znn and Cm transitions we conclude that the depletion in 
these absorbers is uniform at the 0.3 dex level, which corresponds 
to deviations < 3<r in relative abundances at the SNR of the spec- 
tra. It is, therefore, unlikely that these absorbers harbour molecular 
hydrogen, as when H 2 is present the depletion is concentrated in 
one or two velocity components. The absorbers studied here have 
low Wg 934 (< 0.7 A), thus higher W 3934 absorb ers may still trace H 2 . 
These results are consistent with Nestor et aD|2008), whose photo- 
ionisation modeling implied that only high W 3934 Can absorbers 
should be associated with molecular hydrogen. If echelle observa- 
tions of higher systems fail to show evidence for molecular 
hydrogen then it may be necessary to reev aluate the assumptions 
behind the modeling of Ne stor et alj J20081) . such as the strength of 
the ionising UV background and the assumption of solar abundance 
ratios. 

(iii) J0517-4410 stands out from the remainder of the absorber 
sample, exhibiting variation in dust depletion which may map to 
the detection of molecular hydrogen or Ci in this absorber. The 
disadvantage of the AOD method is that binning in velocity space 
smears out the signal of dust depletion, precluding identification of 
the velocity components responsible. Thus it may be that the nar- 
row C 1 and H 2 lines occur at the same velocities, whilst the broader 
Zn n and Cr 11 lines do not correspond to the same gas. Narrow com- 
ponents of Zn 11 or Cr 11 would be hidden by the broader absorption 
profile, so even a detailed Voigt profile analysis of this absorber 
would not illuminate the situation further. Furthermore, ionization 
effects on the meas ured variation in d ust depletion across the profile 
may be important dOuast et al.ll2008h . 

(iv) There are only a few systems where we can study the chem- 
ical uniformity of the absorbers by examining their [Fe/Si] ratios. 
Jl 129+0204 shows ~ 0.3 dex (3<r significance) variation in its 
[Fe/Si] profile that is unlikely be explained by variable dust de- 
pletion, though at the SNR of the data this cannot be certain. 

(v) The Can absorbers studied here are not comprised of ex- 
tremely dense material. Assuming a Milky Way-like radiation field 
there electron densities are, n e < 0.1cm -3 . Thus, unless the gas 
is very neutral (n H +/«H i$ 0.01), n H ~ latomcnT 3 . Comparing 
the three absorber s which overlap between our sample and that of 
Nest or et al. we find a range of UV radiation field strengths 
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are likely in these Can absorbers, from 0.5x to 2.6x the mean 
Galactic field strength. 

(vi) We conclude that most, if not all, Can absorbers studied 
here are not simple disc-like systems because they have complex, 
broad velocity profiles. Simulations only partially constrain the as- 
trophysics of the observed velocity profiles of absorbers as multiple 
mechanisms can produce similar velocity profiles. 

Whilst this work offers the first opportunity to examine Can 
absorbers at high resolution, many questions as to the origin of 
these rare absorbers and the environments of their host galaxies 
remain unanswered. For further ground to be made in our under- 
standing of Can absorption systems, future studies will require 
echelle observations of larger samples, including the, as yet un- 
probed, higher W^ 934 systems as well as host galaxy Integral Field 
Unit (IFU) observations. 
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APPENDIX A: ABSORBER VOIGT PROFILE FITS 

Here we present the Voigt profile fits to our absorption profiles for 
the absorbers and species listed in Table[3] 



APPENDIX B: UNFITTED DATA 

Here we present the Mg n absorption line profiles which were not 
fitted (See Fig. lBlt . We also collate the observed Can profiles into 
one figure as a comparative aid (See Fig. lB2b . 

This paper has been typeset from a TgX/L3TjX file prepared by the author. 
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Figure Al. Voigt profile fits to the z = 0.59760 absorber towards J0334— 0711. A label (bottom-left) in each panel specifies the transitions incorporated in 
the profile fits to the spectral segment shown. The data are represented by the black histograms, whilst the continuum is the dashed line. The fit is shown as a 
solid line through the data. The residuals between the spectra and the fits, normalised by the lcr error, are plotted above the spectra. Individual Voigt profile 
components are marked with ticks. In regions with blends the higher wavelength transition is marked with dashed ticks. The region representing At'90 is lightly 
shaded behind the data. 
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Figure A3. Voigt profile fits to the z = 1.02316 absorber towards J0953+0801, see Fig. lAll for description. 
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Figure A4. Voigt profile fits to the z = 0.83460 absorber towards J1005+1 157, see Fig. lAll for description. 
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Figure A5. Voigt profile fits to the z = 0.74030 absorber towards Jl 107+0048, see Fig. lAll for description. 
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Figure A6. Voigt profile fits to the z = 0.96497 absorber towards Jl 129+0204, see Fig. lAll for description. 
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Figure A7. Voigt profile fits to the z = 0.74630 absorber towards J1203+1028, see Fig. lAll for description. Note there is a blend from an absorber at z. = 0.739 
in the Fen /12374 profile and blends from absorbers at z = 1.32 and 1.57 in the Cm and Znii/12062 profiles. 











1 1 1 1 1 1 1 1 II 1 1 1 1 1 1 




l2 


/-A 


1^ 



Znll2026 Mql2026 



E 0.9 
0.8 













Aft 


Crll2062 


mm 1 1: 1 ml: in: : 

Znll2062 * 





Velocity [km/s] 



Velocity [km/sj 



Velocity [km/s] 




,.2 ^Mfl/vw|^^ 1.2 ^A A ^^ 




Till3242 



-tuu — IUU U , IUU £ 

Velocity [km/s] 



1.0 
0.8 E- 
0.6 E- 
0.4 
0.2 
0.0 




SNH808 



_ U „ IUU <£uu - cuu - IUU , U , IUU 

Velocity [km/s] Velocity [km/s] 




-200 



Velocity [km/s] Velocity [km/s] Velocity [km/s] 

Figure A8. Voigt profile fits to the z = 1.24180 absorber towards J1430+0149, see Fig. lAll for description. 
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Figure A9. Voigt profile fits to the z = 1.14955 absorber towards J0517— 4410, see Fig. lAll for description. 



22 B. J. Zych et al. 




J1430+0149 



-200-100 100 200-100 100 200-200-100 100 200-100 100 200 

2 

°l 

8 

6 



-250-1000 100 250 



J1203+1028 



n 



0.4 



fh*o 



LLLL 



J1 13007- 144927 



■100 100 250-250 -100 100 250 

i.2 r 



100 100 250 




-250 -100 100 250 



100 100 250-250 -100 100 250 
Relotive Velocity fkm/sl 



100 100 250 



-250-100 100 250 



100 100 250-250 -100 100 250 
Relotive Velocity fkm/sl 



100 100 250 




-300-150 150 300-150 150 300-300-150 150 300-150 150 300 
Relotive Velocity [km/s] 

Figure Bl. The Mg n absorption line structure in each of our absorbers with UVES or HIRES data with sufficient spectral coverage. The absorber's quasar 
is identified in the top left-hand corner of each panel. The black histograms are the data and the red histograms are the associated error, whilst the adopted 
continuum is shown by the dashed line. 
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Figure B2. The Ca n absorption line structure in each of our absorbers with UVES or HIRES data with sufficient spectral coverage. The absorber's quasar 
is identified in the top left-hand corner of each panel. The black histograms are the data and the red histograms are the associated error, whilst the adopted 
continuum is shown by the dashed line. 



